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Casimir Force between Atomically Thin Gold Films 
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Abstract. We have used density functional theory to calculate the anisotropic dielectric functions for 
ultrathin gold sheets (composed of 1, 3, 6, and 15 atomic layers). Such films are important components in 
nano-electromechanical systems. When using correct dielectric functions rather than bulk gold dielectric 
functions we predict an enhanced attractive Casimir-Lifshitz force (at most around 20%) between two 
atomically thin gold sheets. For thicker sheets the dielectric properties and the corresponding Casimir 
forces approach those of gold half-spaces. The magnitude of the corrections that we predict should, within 
the today's level of accuracy in Casimir force measurements, be clearly detectable. 
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1 Introduction 

Casimir pQ predicted already in 1948 that boundary ef- 
fects on the electromagnetic fluctuations can produce at- 
traction between a pair of parallel, closely spaced, per- 
fectly conducting plates. This work was later extended 
to real materials by Lifshitz [2J|3] , and by Parsegian and 
Ninham \4\ . Since the famous experiments of Deryaguin 
and Abrikossova [5] there has been much interest in phe- 
nomena that measure the van der Waals forces acting 
between macroscopic bodies. The early experiments that 
measured the forces between quartz and metal plates cov- 
ered only the retarded region. The experiments of Ta- 
bor and Winterton [5] and subsequently of Israelachvili 
and Tabor [7J|8] fitted the potential to a power law of 
l/d n (d being the distance) where n varied from non re- 
tarded (n = 3) to fully retarded (n = 4) value. There 
was a gradual transition from non-retarded to retarded 
forces as the separation was increased from 120 A to 1300 
A [7] . Lamoreaux performed the first high accuracy mea- 
surement [§] of Casimir forces between metal surfaces in 
vacuum [31110] . The first measurements of Casimir-Lifshitz 
forces applied to micro-electromechanical systems were 
performed by Chan et al. [T111T2] and somewhat later by 
Decca et al. [13J . Several new very recent high precision 
measurements have been performed fiTLIMTSlITT] . An in- 
teresting aspect of the Casimir-Lifshitz force is that ac- 
cording to theory it can be either attractive or repulsive [21 
fT51fl9 . 20,21,22. 23 . Casimir-Lifshitz repulsion was mea- 
sured four decades ago for films of liquid helium (10-200 
A) on smooth surfaces [24] . Only a few direct force mea- 
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surements of repulsive Casimir-Lifshitz forces have been 
reported in the literature [251126112711251125] . 

Surfaces that are very close interact in vacuum with 
a van der Waals force [Tlf^f5lf6lf30l[3T] [32 l. f33 | f34 |f35l[36] . The 
van der Waals force between thin isotropic metal films fol- 
lows a fractional power law related to the two dimensional 
plasmon dispersion [37U38U39U4D1I41U18] . As the separation 
between the surfaces increases the interaction takes a weaker 
(retarded) Casimir form [TJ|3] . As demonstrated by Be- 
nassi and Calandra [12] the Casimir force between ultra- 
thin (1 to 10 nm thick) conducting films depends on the 
anisotropy of the films. In the recent Casimir experiments 
the metal films coating the objects were thick enough 
for the objects to be considered as bulk metal objects. 
The metal films were assumed to be homogeneous and 
isotropic. This is not always so in reality. Svetovoy et 
al. [43] measured both the optical properties and Casimir 
forces for 1000-4000 A thick gold films on substrates. As 
a reference curve for the dielectric function of thin films 
they |45] chose one which was calculated with handbook 
data |33] . The Casimir force evaluated for their films was 
considerably weaker compared to that calculated with the 
reference curve (5% to 14% weaker at a distance of 100 
nm between the films). Thus, the dielectric function and 
force are sensitive to how the films are prepared. Thin 
films of gold can even be insulating [45] due to disorder. In 
Reference [31] one calculated the force near the insulator- 
conductor transition in thin gold films. In the present work 
we are not concerned with these effects. We assume homo- 
geneous ordered metallic films. We address what happens 
when the metal films are extremely thin. 

We first present the theory used to calculate anisotropic 
dielectric functions of atomically thin gold films. Such 
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films are important components in nano-electromechanical 
systems (NEMS). Then we demonstrate that there is an 
enhancement of the Casimir-Lifshitz force (of the order 
of 20 %) when proper dielectric functions for ultrathin 
films are used rather than the dielectric function appro- 
priate for bulk gold. As the film thickness increases the 
use of isotropic bulk dielectric function becomes an in- 
creasingly good approximation. The long range entropic 
Casimir asymptote originates entirely from the transverse 
magnetic zero frequency mode in agreement with recent 
experiments |15j . We end with a brief summary and an 
outlook. 



2 Calculation of dielectric permittivities of 
atomically thin gold films 

We have calculated the anisotropic dielectric function of 
ultrathin gold sheets (1, 3, 6, and 15 atomic layers) plus 
the corresponding results for thick gold plates. The calcu- 
lations of the dielectric function were performed within the 
density functional theory, employing the local density ap- 
proximation in conjunction with the projector augmented 
wave method 17. 1?*. H7]. We modeled the layer structures 
by supercells oriented in the crystalline (111) direction. 
The supercells have hexagonal symmetry with a height of 
56.5 A(i.e., 24 atomic layers). The imaginary part of the 
dielectric function was calculated in the long wave length 
limit from the optical momentum matrix, and the intra- 
band contribution was included by the Drude-like model 
with a damping of 7 = 5 meV. The dielectric function 
on the imaginary frequency axis was determined from the 
Kramcrs-Kronig dispersion relation. The ratio of these di- 
electric functions to the corresponding dielectric function 
of bulk gold are shown in Figure[T] (where z is the (111) 
direction perpendicular to the thin film and x is a direc- 
tion within the film). These dielectric functions have been 
used together with the expression for the anisotropic Lif- 
shitz force given by Benassi and Calandra [42] , with proper 
extension to include finite temperature effects [3D] ■ While 
Casimir forces between thin films (e.g. lipid films in water) 
have been known for more than 40 years [30 ,35,50 they 
have not been explored for the present case of ultrathin 
anisotropic conducting films. Nanotechnological advances 
and refined atomic models of ultathin gold sheets now al- 
low their exploitation. 

3 Theory and Numerical Results 

The simlest way to find retarded van der Waals or Casimir- 
Lifshitz force is in terms of the electromagnetic normal 
modes 30,32, 42 of the system. At finite temperature the 
force (F) originates from changes in the Helmholtz' free 
energy and can be written as 

OO r o 

J ° „^ !K UJ ll-QTMi^n) 2 (1) 
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Fig. 1. (Color online) Ratio between the diagonal elements 
of the dielectric tensor of ultrathin gold sheets and that of 
bulk gold. The in-plane elements, e xx = e yy , (a), are enhanced 
compared to the bulk value while the perpendicular component 
£zz, (b), is reduced in value. e zz (iui) = 6.4, 21.9, 60.2, and 
198 for sheets composed of N — 1, 3, 6, and 15 atomic layers, 
respectively, whereas bulk gold has a value of e(iuo\) = 2700. 
The lowest nonzero frequency ui\ is ~ 2.47 x 10 14 rad/s. 
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Fig. 2. (Color online) Casimir-Lifshitz force between ultrathin 
gold films, between gold half-spaces, and between ideal (per- 
fectly reflecting) metal surfaces. For comparison we also show 
the n = entropic long range asymptotic force, which is the 
same for all cases except for the Casimir ideal case. 
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(2) 

(3) 

(4) 
(5) 
(6) 

(7) 

In Equation ([2]) I is the distance between the surfaces 
of the two sheets and D = [2e?o + (N — l)d] is the film 
thickness ( do = d/2, d = 2. 35 A and N is the number of 
atomic layers) . With this definition the thickness of a sin- 
gle layer is the same as the distance between two atomic 
layers. The integral over frequency appropriate for zero 
temperature [42) has been replaced by a summation over 
discrete Matsubara frequencies to account for finite tem- 
perature. The prime on the summation sign indicates that 
the n = term should be divided by two. For planar struc- 
tures the quantum number that characterizes the normal 
modes is k, the wave vector in the plane of the interfaces, 
and there are two mode types, transverse magnetic (TM) 
and transverse electric (TE). 

We show in Figure[2]the Casimir-Lifshitz force between 
atomically thin gold sheets (exemplified with the cases of 
1 layer and 15 layers). All results are for T = 300K. The 
result is compared with the Casimir-Lifshitz force between 
two gold half-spaces and between two ideal metal sur- 
faces. Included is also the long-range entropic contribution 
(n = term in the Matsubara summation). This is the 
same for all cases except for the Casimir ideal case where 
this contribution is twice as big. For the range covered 
in the figure temperature effects have not yet become im- 
portant. The force between two gold half-spaces is weaker 
than the force between two ideal metal half-spaces. It ap- 
proaches this force at the rightmost end of the figure where 
retardation effects prevent the crossing of the curves. The 
force between two 15-layer films approaches the force be- 
tween two gold half-spaces well before the leftmost end of 
the figure is reached. For the 1-layer films this happens 
much later, outside the figure. 

In order to study the effect of using dielectric func- 
tions calculated for thin films we present in Figure[3J the 
ratio between Casimir-Lifshitz force calculated with real- 
istic anisotropic dielectric functions to the corresponding 
calculation using isotropic bulk dielectric function. We ob- 
serve corrections in the range between 5 and 20 % in the 
range most accessible to force measurements. 
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Fig. 3. (Color online) Ratio of the Casimir-Lifshitz force be- 
tween thin films calculated with realistic anisotropic dielectric 
functions and from the corresponding force calculated using 
the isotropic bulk dielectric function. 
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Fig. 4. (Color online) Ratio of Casimir-Lifshitz force between 
thin films and corresponding force between gold half-spaces. 



We finally demonstrate in Figure|3J how the Casimir- 
Lifshitz force between thin gold films is reduced compared 
to the corresponding force between gold half-spaces. 

To be noted is that the maximum deviation from unity 
in both Figures[3]and[3]occurs in the middle portion of the 
figures. For large separations only small frequencies con- 
tribute. For zero frequency Ptm = ~ 1 and the factors 
containing the film thickness in Equation @ cancel. Fur- 
thermore pte = and Qte — for zero frequency. This 
behavior results from the properties of the low frequency 
screening: e xx (uS) goes towards infinity at low frequencies 
and u) 2 e xx (ui) goes towards zero. For small separations the 
force between two films approaches the force between two 
half-spaces [55] and the anisotropy effects are reduced. 
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4 Conclusions 

It has been stressed how crucial it is in Casimir force 
calculations to use accurate dielectric functions from opti- 
cal data or from calculations. We demonstrate in this arti- 
cle that when performing calculations on Casimir-Lifshitz 
forces between ultrathin conducting films it is not suffi- 
cient to use dielectric functions relevant for thick films or 
half-spaces. Rather we have demonstrated the importance 
of using proper anisotropic dielectric functions from den- 
sity functional calculations when performing calculations 
of Casimir forces between atomically thin gold sheets. 
These ultrathin films have interesting applications in nano- 
technology, including use in NEMS. Addition of ultrathin 
metallic coatings to surfaces in solution have for instance 
been predicted to create quantum levitation of NEMS [22 . 
It remains to be seen how such predictions may be influ- 
enced when using accurate dielectric functions of atomi- 
cally thin metal films. 

CP. and M. B. acknowledge support from VR (Contract No. 
90499401 ) and STEM (Contract No. 34138-1). B.E.S. acknowl- 
edges financial support from VR (Contract No. 70529001). 
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